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Abstract 

O . 

This is a didactic review of spectral and dynamical properties of atoms and molecules 
£SJ ■ at energies below the ionization threshold, the focus being on recent work in which the 

author was involved. As far as possible, the results are described using a simple model 
with one electron only, and with scalar bosons. The main ideas are explained but no 
complete proofs are given. The full-fledged standard model of non-relativistic QED and 
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C*~) , various of its aspects are described in the appendix. 
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An atom or molecule in an excited state with energy below the ionization threshold will eventu- 
ally relax to its ground state by dissipating excess energy in the form of radiation. This process 



43 , 

of relaxation to the ground state is one of the basic phenomena responsible for the production 
of all visible light. It involves a range of energies within a few electron volts; a scale where the 
electron-positron pair creation and the production of ultraviolet radiation is highly suppressed. 
In a first mathematical study of relaxation to the ground state it is therefore reasonable and 
legitimate to work with a model where the electron-positron pair creation is entirely neglected 
and an ultraviolet cutoff is imposed on the electron-photon interaction. These simplifying as- 
sumptions lead to a mathematically well defined model of matter, often called standard-model 
of non-relativistic quantum electrodynamics, or Pauli-Fierz model. Since the numerical predic- 
tions of this model are in good agreement with measurement data it is a viable physical model. 
Yet, only little mathematically rigorous work on this model had been done before the middle 
of the 1990s, when several groups of researchers started to investigate various of its aspects. 
Most influential, perhaps, were the papers of Hiibner and Spohn [331 dH 1351 on spectral and 
scattering theory, of Bach et al. on spectral analysis [71 1§1 151 HTJj . of Dereziiiski and Gerard on 
scattering theory [141 115j , and of Jaksic and Pillet on thermal relaxation [371 1381 1391 14U| . The 
present article reviews recent work on the phenomenon of relaxation to the ground state for 
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states with total energy below the threshold energy for ionization. It is guided by papers of 
Lieb, Loss and Griesemer, and of Frohlich, Schlein and Griesemer |251 1241 IT71 118j . The focus 
is on the existence of an ionization threshold and the localization of the electrons below this 
energy, the existence of a ground state, and the existence and completeness of many-photon 
scattering states (asymptotic completeness of Rayleigh scattering). 

The results to be discussed on existence of a ground state and on the localization of pho- 
tons with energy below the ionization threshold, unlike previous results, hold for all values of 
the physical parameters such as the fine structure constant and the ultraviolet cutoff. This 
is crucial for moving on to physically more realistic models without ultraviolet regularization. 
The analysis of electron-photon scattering is based on methods and ideas from the scattering 
theory of iV-body quantum systems |5U1 1231 1541 122 j . On the one hand the electron-photon 
dynamics is easier to analyze than the full A-body problem since there is no photon-photon 
interaction. On the other hand the number of photons is not constant! In fact, it might even 
diverge as time t — ► oo. This divergence is avoided by imposing a cutoff on the interaction 
between electrons and low-energy photons (infrared-cutoff). It is one of the main open chal- 
lenges in the mathematical analysis of matter interacting with quantized radiation to prove 
asymptotic completeness for Rayleigh scattering without this infrared-cutoff. 

This article is organized as follows. In Section 2.1 we begin with the description of a simple, 
but non-trivial model of matter and radiation. There is only one electron, besides the static 
nuclei, and the radiation is described by scalar bosons. 

Sections 12.41 12.51 12.71 and 12.81 describe the main results of the papers j^J HSJ El and 
|18| . respectively. Sections 12 . 21 and 12 . 31 summarize mathematical and physical background, and 
Sections 12.91 and 12.101 are devoted to side issues in the aforementioned papers. 

Section outlines the modification of results and proofs that are necessary to accommodate 
N > 1 electrons, and Section 0] ends this review with concluding remarks and a discussion of 
selected open problems. 

There is a self-contained appendix on the standard model of non-relativistic QED. 

Acknowledgments. Most of the content of this review article I learned from my collaborators 
Jiirg Frohlich, Elliott H. Lieb, Michael Loss, and Benjamin Schlein. I am indebted to all of 
them. I thank David Hasler for his careful proofreading. 

2 Matter and Radiation 

All electrons of an atom or molecule are well localized near the nuclei if the total energy is below 
the ionization threshold. Therefore the number of electrons is inessential for the phenomena 
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to be described mathematically in this section. To simplify notation and presentation we 
restrict ourselves to one-electron systems; the generalization to N > 1 electrons is described 
in Section |21 

2.1 A Simple Mathematical Model 

The main features of quantum electrodynamics that are responsible for the phenomena to be 
studied, are the peculiar form of interaction between light and matter, through creation and 
annihilation of photons, and the fact that photons are massless relativistic particles. The spin 
of the electron and the helicity of the photons do not play an essential role in most of our 
analysis. For the purpose of this introduction we therefore neglect these subtleties and present 
a caricature of QED which only retains the aforementioned main features. The full-fledged 
standard model is described in the appendix. 

We first introduce our models for matter and radiation separately before describing the 
composed system and the interaction. 

A (pure) state of a quantum particle, henceforth called electron, is described by a normalized 
vector ip G L 2 (R 3 ,C), J, \ip(x)\ 2 dx being the probability to find the particle in the region 
AcR 3 . Its time evolution is generated by a Schrddinger operator 

flat = -A + V (1) 

where —A is the positive Laplacian and V is the operator of multiplication with a real-valued 
function V(x), x G M 3 . We assume that V G L 2 oc (M 3 ) and that there exist constants a < 1 
and such that 

{<p,V-<p)<a{<p,i-A)<p)+P{<p,tp) (2) 

for all ip G Cq°(R 3 ), where V— = max(— V, 0). Hence the operator —A + V is symmetric and 
bounded from below, which allows us to define a self-adjoint Hamiltonian fl at by the Friedrichs' 
extension of — A + V. 

We shall be most interested in the case where 

V(x) = V z (x):=-J2 ] -§—, (3) 

j=l 3 1 

Zj, j = 1 . . . K, are positive integers, and Rj G M 3 . The function Q is the potential energy 
(or the scalar potential in Coulomb gauge) of one electron at x G M 3 in the field of K nuclei 
with positions R\, . . . , Rk and atomic numbers Z\, . . . , Zk- 

The Hamiltonian (JJ) with V given by © describes a molecule with one electron and 
static nuclei in units where the unit of length is h 2 / (2me 2 ) = re/2 and the unit of energy is 
2e 2 /re = 4 Ry (see Appendix lA")) , Here r# = H 2 /(me 2 ) is the Bohr radius, — e is the charge of 
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the electron and m is its mass. This Hamiltonian is self-adjoint with domain D{H &t ) = ff 2 (lR 3 ), 
the Sobolev space of twice weakly differentiable L 2 -functions [1T1 H7j . 

A pure state of the radiation field is described by a normalized vector in the bosonic Fock 
space over L 2 (M 3 ). This is the space 

f=05 n L 2 (l 3n ;C) 

n>0 

where SoL 2 (M°) := C, and S n denotes the orthogonal projection onto the subspace of square 
integrable functions f(k\, . . . ,k n ) that are symmetric with respect to permutations of the n 
arguments k%, . . . , k n G M 3 . Such a function describes a state of n bosons, henceforth called 
photons, with wave vectors ki, . . . , k n . The vector |vac) = (1, 0, 0, . . .) G T is called the vacuum 
vector. With J 7 ^ we denote the subspace of sequences cp = (ip) n >o G J 7 with <p n = for all 
but finitely many n G N. 

The energy of a state (p = (ip n )^ =0 G T is measured by the Hamiltonian Hf defined by 

(H f <p) = 

" (4) 

(H f ip) n (ki, ...,k n )= 2_^u{kj)<Pn(kl, ■ ■ ■A), Tl > 1, 

i=i 

where uj(k) = \k\. The domain of Hf is the largest set of vectors for which Q) defines a vector 
in J-. 

The interaction between photons and electrons comes about in a process of creation and 
annihilation of photons. To describe it mathematically, creation and annihilation operators 
are needed. Given h G L 2 (M 3 ) and <p G !F^ a we define a*(h)ip by 

[a*(h)<p] n = ^S n (h (g> <p n -i)- 

The operator a*(h) is called a creation operator. It adds a photon with wave function h to the 
state The annihilation operator a(h) is the adjoint of the closure of a*(h). These operators 
satisfy the canonical commutation relations 

[a(g),a*(h)] = (g,h), [J(g), a*(h)] = 0. (5) 

A further important operator on T is the number operator Nf, defined by 

(Nf<p) n = rup n 

and D(Nf) = {<p G T : E n2 |^n|| 2 } < oo. 

A state of the composed system of electron and photons is described by a vector ^ G 
Wat <S5 T ', that is, by a sequence (ip n )^ =0 where ip n is a square integrable function 

ip n (x, k±, . . . , k n ), 
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describing a state of one electron and n photons. It is often helpful to use that Ti. a t <8> T — 
L 2 (M 3 ;^ r ) and to consider 'J' as a square integrable function x i— > ^(x) with values in T. Then 
||^ r (x)||^r is the probability density for finding the electron at position i£l 3 , 
For the generator of the time-evolution t i— ► *$> t we choose the Hamiltonian 



where j£l and k € Co°(M 3 ). It is easy to prove that H[ nt is operator-bounded with respect 
to H g= Q with bound zero, and hence, for every g € R, i? s is bounded below and self-adjoint 
on the domain of H g= o by the Kato-Rellich theorem |47j . 

Both (7 and k measure the strength of interaction between electron and photons. Given k, 
some of the results in the following sections hold for \g\ small enough only, others for small 
\g\ 7^ 0. The value n(k) of the form-factor k measures the strength of interaction between 
electron and radiation with wave vector k. There is no interaction for k outside the support 
of k, which is the case, e.g., for \k\ larger than the ultraviolet cutoff A := sup{|/c| : n{k) 7^ 0}. 

We reiterate that this model is a caricature of the standard model of quantum electrody- 
namics for atoms and molecules interacting with quantized radiation. With QED it has in 
common that it describes a non-relativistic particle, the electron, interacting with massless 
relativistic bosons in a momentum conserving process of creation and annihilation of such 
bosons. 

It would make our toy model physically more realistic if we assumed K,(k) ~ |/c| -1 / 2 for 
small \k\. But then H has no ground state [11], a problem that does not occur in the standard 
model of non-relativistic QED. Therefore we assume that k is non-singular near k = 0. 

2.2 Spectrum and Eigenfunctions of H at 

As a preparation for the following sections we recall a few facts concerning the spectrum of 
Schrodinger operators H &t = —A + V and the decay of their eigenfunctions. 

Suppose that V satisfies assumption (J2J) and let H^t be defined in terms of the Friedrichs' 
extension of the symmetric operator —A + V on C^°(M 3 ). Let Dr = Co°(|x| > R), the space 
of smooth, compactly supported functions with support outside the ball Br(0), and let 



H = H g = H at (8) 1 + 1 (8) H f + gH int 



(6) 



the interaction H mt being given by 



(H- mt V)(x) = [a(G x ) + a*(G x )Mx) 
G x (k) = e-^ x K(k), 




(7) 
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By a theorem due to Arne Persson |46l [2] 

S at = inf cr ess (iJ at ) (8) 

where a ess (H a t) denote the essential spectrum of H at , i.e. the complement, within the spec- 
trum, of the isolated eigenvalues of finite multiplicity. For the Coulomb potential ©, Vz{x) — ► 
as |x| — > oo and hence S a t = 0. Furthermore 0" ess (i7 a t) = [0, oo), and by a simple variational 
argument, -ff a t has infinitely many eigenvalues below |48j . 

Eigenfunctions of H^t with energy below S a t decay exponentially with increasing \x\: for 
every eigenvalue E < S at and every (3 > with E + (5 2 < S at there exists a constant Cp such 
that 

\ij}(x)\ < Cpe-PM, a.e. on M 3 (9) 

for all normalized eigenfunctions that belong to E. Of course, the actual decay of tp will not 
be isotropic unless V is spherically symmetric. There is the better, but non-explicit, bound 
|i/>(x)| < C , e e - ( 1-£ '' p ( :c ) where p(x) is the geodesic distance from x to the origin with respect to 
a certain metric ds 2 = ce(x/\x\)(1x 2 in M 3 Since ce{x/\x\) > S at — E (if S at < oo), the 
isotropic bound © follows from this stronger result. 
For proving © it suffices to show that 

e^ 1 ?/' G L 2 (R 3 ) (10) 

whenever E + (3 2 < inf a ess {H a \). The point- wise bound Q then follows from a general result 
on point-wise bounds for (weak) solutions of second order elliptic equations [21 The L 2 - 
bound (|l(Jj) . in turn, is easily derived from the characterization (J7J), © for inf a ess (H a t) 

We now turn again to H g , the Hamiltonian @ describing matter and radiation. Most 
properties of H g to be discussed in the following sections hold for all j £ 1, including g = 0. 
Hence they generalize properties of 

H g=0 = H at ® 1 + 1® H f . 

By general spectral theory o{Hq) = a(H at ) + a(Hf). Furthermore, it is easy to see from the 
definition of Hf that cr(Hf) = [0, oo) and that is the only eigenvalue of Hf, the vacuum |vac) 
being its eigenvector. It follows that cj{Hq) = [inf c(if a t), oo) and that Hq and H a t have the 
same eigenvalues. The corresponding eigenvectors are the products tp (g) |vac), where ip is an 
eigenvector of H at . 

2.3 Physical Phenomena and Mathematical Description 

The experimental evidence on isolated atoms in contact with radiation is easiest described in 
an idealized setup or "Gedankenexperiment" . Consider an atom in a universe that is otherwise 
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free of matter. For simplicity, we assume that the atom has only one electron. There may 
be radiation near the atom initially but no "external fields" or sources of radiation shall be 
present. Independent of its initial state, this system of atom and radiation will eventually 
approach one of only two qualitatively distinct final states: the "bound state" or the "ionized 
state" . 

In the bound state the atom is in its ground state, the state of least energy, where the 
electron is confined to within a small neighborhood of the nucleus. All excess energy has been 
radiated off. This radiation is very far away from the atom and escaping at the speed of light. 

In the ionized state the electron and the nucleus are spatially separated with increasing 
distance. In addition, there may be radiation going off to infinity. 

Of course, the ionized state can only be attained if the total energy of matter and radia- 
tion initially is hight enough to overcome the attraction between nucleus and electron. High 
energy however, does not guarantee ionization, as the excess energy may just as well turn into 
radiation. On the other hand, if the total energy initially is not sufficient for ionization, then 
the atom will certainly relax to its ground state. Other conceivable scenarios, like relaxation 
to a stationary state with non-minimal energy, or the permanent radiation without total loss 
of the excess energy have never been observed. Bohr's stationary states, with the exception 
of the ground state, are unstable, and radiation is only emitted in the very short period of 
transition to the ground state. 

The goal is to give a proof of the phenomenon of relaxation to the ground state in the model 
introduced in the previous section. In view of the experimental evidence described above, we 
expect that this model has the following mathematical properties. 

Existence of the ionization threshold. There exists a threshold energy S > 'mta(H), 
such that the electrons described by states in the spectral subspace -E7(_oo,£)(-£0^ are 
well localized near the nuclei, while states with energy above S may be ionized. 

Existence of a ground state. There exists a state of least energy (ground state), or, equiv- 
alently, inf o~(H) is an eigenvalue of H . 

Absence of excited stationary states. The operator H has no eigenvalues above inf cr(H). 

Asymptotic completeness of Rayleigh scattering (ACR). In the limit t — ► oo the time 
evolution e -tfft \E r of every state ^ 6 Ei 00 ^){H)'H is well approximated by a superposi- 
tion of states of the form 



where ha(k) = e ^ k ^hi(k), is a normalized ground state of H and Eq is its energy. 
The vector (flip describes a state composed of the atom in its ground state and n freely 
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propagating photons with wave functions hi t, . . . , h n In the limit t — > oo they will be 
far away from the atom. 

The papers HUHSIEIEI are devoted to proving the above four properties for the standard 
model of QED was well as for the model introduced in Section fe.ll One exception is asymptotic 
completeness for Rayleigh scattering where the results concern the model of Section [2.11 only, 
with the additional important simplification that n(k) = for \k\ < a where a > is arbitrarily 
small but positive. This assumptions that photons with energy near zero don't interact with 
the electron is usually referred to as an infrared cutoff, IR-cutoff, for short. Sometimes the 
constant a is called infrared cutoff as well. The significance of the IR-cutoff is that it allows 
us to control the number of bosons that are being produced in the course of time: the number 
of bosons with energy below a stays constant and the number of bosons with energy above a 
can be bounded from above in terms of the total energy. The assumption of an IR-cutoff is not 
expected to be necessary for the validity of asymptotic completeness of Rayleigh scattering as 
formulated above; but as of now, no convincing mathematical argument is known that would 
substantiate this believe. 

The above list of physical phenomena is limited to the coarsest properties of atoms inter- 
acting with radiation. Even below the ionization threshold, there are many other phenomena 
that are worth rediscovering in the standard model, and in part this has already been done. 
Most important, perhaps, are the occurrence of sharp lines in the spectrum of the emitted ra- 
diation (Bohr frequencies), the resonances in lieu of Bohr's stationary state and their extended 
life time, and the correspondence principle at energies near the ionization threshold. We shall 
come back to some of these phenomena in later sections, when we review known results or 
comment on open problems. One must keep in mind, however, the limitations of our model. 
Quantitative predictions will be of limited accuracy when relativity or high-energy photons 
play a significant role. 

2.4 Exponential Decay and Ionization Threshold 

The ionization threshold £ of an atom or molecule with only one electron is the least energy 
that this system can achieve in a state where the electron has been moved "infinitely far away" 
from the nuclei. The electron is outside the ball \x\ < R with probability one, if its wave 
function ^f(x) vanishes in this ball. Therefore we define 





where 



D R = {tf e D(H)\V(x) = if \x\ < R}. 
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Note the analogy with Persson's characterization (JJJ) of inf a ess (H at ). Here, however, £ ^ 
inf a ess {H ) unless S = inf cr(H). In general £ > inf cr(-ff) and, if V(x) —* as |x| — > oo, then 
£ = inf <r(-ff — V), which is greater than inf cr(H) for V = Vz ( see also Section [2. 5p , 

According to (|12|) . the electron described by a state \& with energy below S cannot be 
arbitrarily far away from the origin. In fact, by [21], for all A, f3 E M with A + f3 2 < S 

|| (e^l ®1)£ A (#)|| < oo- (13) 

This shows that the probability(-density) ||^(x)|||r to find the electron at the point x decays 
exponentially fast, as |x| — » oo, at least in the averaged sense 

j e W \\^{x)\\%dx < oo. (14) 

There is an obvious similarity between 1)141) and 1)10(1 that is not accidental. In |24| the 
bound 1)13)1 is derived from an abstract result for semi-bounded self-adjoint operators H in 
Hilbert spaces of the form L 2 (IR n ) <g> where T is an arbitrary complex Hilbert space. The 
only assumptions are that fD(\H\ 1 / 2 ) C D(\H\ 1 / 2 ) and that 

f 2 H + Hf 2 -2fHf = -2\Vf\ 2 (15) 

for smooth, bounded functions f(x) with bounded first derivatives. Equation (|15|) holds for 
H = — A <S> 1 and since the left-hand side of ()15|) formally equals [/, [/,!?]], it follows that 
()15)) holds for are large class of self-adjoint operators H whose principal symbol is given by the 
Laplacian. The result ifTUJl thus emerges as a special case of p3|). 

Even though the above assumptions on /7 are largely independent of Hf and gH{ nt , the 
result depends on these operators! The binding energy S — inf cr(H) depends of g and hence 
so does the decay rate one obtains for the ground state. 

It is well known, since the work of Agmon |2], that eigenfunctions of second order elliptic 
equations decay exponentially in energetically forbidden regions, that is, in regions where the 
differential operator, as a quadratic form, is strictly larger than the eigenvalue [2]. The result 
()13)) shows that this idea can be brought to bear in a much more general framework, including 
many models of non-relativistic QED. It is clear that our proof of ()13|) can be generalized, along 
the lines of |2], to yield non-isotropic bounds, as well as exponential bounds for QED-models 
where H^t is a more general, uniformly elliptic second order differential operator. 

We owe the strategy for proving ()13|) to Bach et al. [7], where this bound is established for 
\g\ sufficiently small and A + f3 2 < S at — const \g\. 

2.5 Existence of a Ground State 

By the main result of j^S], infer (H) is an eigenvalue of H whenever 

inf a(H) < E. (16) 
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This reduces a difficult spectral problem to a variational problem: the problem of finding a 
state 6 D(H) with ii^) < £(\l/, There are two important classes of potentials V for 
which this variational problem can be solved without much effort: if V(x) -> oo as \x\ — ► oo, 
then obviously S = oo and hence inf cr{H) < S. On the other hand, if V(x) — > as |x| — > oo 
then £ = inf a[H— V) where the nuclei are removed in H — V. Using the translation invariance 
of H — V one shows that 

inf a(H) < inf a(H — V) + inf a(-A + V). 

It follows that inf cr(H) < £ whenever V(a;) — > 0, (|x| — ► oo) and inf a(— A + V) < 0. Since 
this is the case for the Coulomb potential V = Vz, all one-electron atoms and molecules have 
a ground state. 

We next sketch the proof that Q16|) guarantees existence of a ground state. To begin 
with we recall that a hypothetical eigenvector ^ of H, with eigenvalue inf cr(H), minimizes 
the quadratic form ^ subject to the constraint ||^|| = 1. It is thus natural to 

establish existence of \E' by proving relative compactness for a suitable minimizing sequence. 
The problem with this approach is that a generic minimizing sequence will tend weakly to 
zero. The fact that inf cr(H) belongs to the essential spectrum alone implies that there are 
infinitely many energy minimizing sequences with this defect. Our task is thus, first, to choose 
a suitable minimizing sequence, and second, to prove its relative compactness. 

We choose the elements of our minimizing sequence to be the ground states \I/ m of modified 



Hamiltonians H m (m — > 0) in which the photon energy u(k) is altered to be oj m {k) = \/k 2 + to 2 . 
That is, we give the photons a positive mass m. For m small enough the binding assumption 
(|TE|) is inherited by H m , which we use to show that inf a(H m ) is indeed an eigenvalue. As 
a matter of fact, inf a{H m ) is separated from the essential spectrum of H m by a gap of size 
m JH]- The sequence of ground states ( l I / m ) m >o is a minimizing sequence for H that can be 
assumed to be weakly convergent. It remains to show that the weak limit is not the zero vector 
in H. 

We first argue that it suffices to prove relative compactness of the sequence of L 2 -functions 
ij} mn {x, k\, . . . , k n ) restricted to large balls B C IR( 3+3n ). This follows from the exponential de- 
cay w.r.to x, from ifi min (x, ki, ■ ■ ■ , k n ) = if \ki\ > A, and from the bound sup m (^ m , Nf^ m ) < 
oo. Then we use the compactness of the embedding 

W l -"(B) •-, L 2 (B), for2>p >!_!±^ 

due to Rellich-Kondrachov p. We thus need to show that sup m ||V^ mjn || p < oo, which we 
derive from a refinement of the argument that we used to prove the bound on ($ m , Nf^f m ), 
and from the -ff-boundedness of — A x . 
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There was a large number of previous papers on the existence of a ground state for models 
similar to the one discussed here @J 13 E3 EM Ell Among these, the best result 
is due to Bach et al 0. It established existence of a ground state when inf cr(if at) jv) is an 
isolated eigenvalue and the fine-structure constant a is small enough. Most importantly, in 
this paper for the first time existence of a ground state is proven in the standard model without 
an IR-regularization. 

2.6 Relaxation to the Ground State is a Scattering Phenomenon 

As discussed in Section 1231 every state ^/ £ RanE^^^^H) is expected to "relax to the ground 
state by emission of photons". In mathematical terms, this means that e~ lHt ^, in the distant 
future, t — > oo, is well approximated in norm by a linear combination of vectors of the form 

a*(h 1>t )...a*(h n , t )e- iEot y (17) 

where hij(k) = e~ lujt hi(k), is the ground state, and Eq is its energy. This has a chance to 
be correct only if H has no other eigenvalues below X. If it has, then relaxation to a bound 
state may occur, which means that i n (|17|) may be any eigenvector of H with eigenvalue 
below S. In this weaker form, the above assertion is called Asymptotic Completeness for 
Rayleigh scattering. The problem of proving absence of excited eigenvalues is independent of 
the scattering problem and its discussion is deferred to a later section. 

Before proving completeness of the scattering states one needs to address the problem of 
their existence. An example of a scattering state is a vector S Ti. for which there exist 
n > 1 photons hi, . . . , h n and an eigenvector ^f, H^f = E*S> such that 

e~ im ^ + ~ a*(h ltt ) . . . a*(h n , t )e- iEt V, t -> oo 

in the sense that norm of the difference vanishes in the limit t — > oo. The scattering state ^+ 
is said to exists if the limit 

* + = lim e iHt a*{h lt )...a*{h nt )e- iEt y t -> oo (18) 

exists. Let TL + denote the closure of the space spanned by vectors of the form (|18|) . All 
elements of TL + are called scattering states and asymptotic completeness of Rayleigh scattering 
is the property that 

W+DRantf^jgGff). (19) 

Existence of scattering states is established in |17j . and (|19f) is proven in |18j for the model 
introduced in Section 12.11 assuming either an infrared cutoff on the interaction or that the 
photon dispersions relation uj(k) is bounded from below by a positive constant, which excludes 
u(k) = \k\. The latter assumption serves the same purpose as the infrared cutoff, and it is 
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satisfied, e.g., for massive bosons where u(k) = \Jk 2 + m 2 with m > 0. As of today, there is no 
proof of ()19|) without a form of infrared cutoff or another drastically simplifying assumption 

CUED- 



2.7 Existence of Scattering States 

Generalizing (|18|). we ask whether the limits 

* + = lim e iHt a#(h ljt ) . . . a*{h n ^ im 



(20) 



exist for given ^ eH, and hi e L 2 (R 3 ), where a*(h i}t ) is a creation or an annihilation operator. 
The scattering states obtained in the limit t — > — oo are physically interesting as well, but 
the problem of their existence is mathematically equivalent to the existence of (|20j) . Beginning 
with the easiest case, n = 1, let us ask whether the limits 



a*{h)t> = lim e iHt a*{h t )e- iHt ^ 



(21) 



t— >oo 



exist. If the photons are massless, as they are in nature, the answer depends on the electron 
dispersion relation and on the energy distribution of ^. For massive photons, however, a^_{h)^ 
exists for all ^ € D(H) and all h G C^°(1R 3 ), and the proof is short and easy |S2] : by the 
Cauchy criterion the limit <\21\i exists if the time derivative of the right-hand side is (absolutely) 
integrable. A straightforward computation using a* (h t ) = e" iH / i a*(/i)e iff / i and © shows that 



where 



^e lHt a*(h t )e~ iHt 
at 



(G x ,h t 



dt 



\(G x ,h t )e- iHt *\\dt 



Ak-x—iuit 



n(k)h(k) dh 



(22) 
(23) 



and uj(k) = Vk 2 + m 2 . Since the Hessian of uj is strictly positive, 

sup \(G X , h t )\< const t~ 3/2 , t > 1 

X 

by a standard result on oscillatory integrals [301 • Hence \\(G X , ht)e~ tHtx &\\ < const t~ 3 / 2 and 
(12'2|) is finite which proves the existence of a*^_(h)^. For massless photons, however, 

sup \(G X , h t )\ ~ const t , t>l 

X 

which is not integrable and we need to estimate the integrand of (|22[) more carefully. To begin 
with, we note that the phase in (J23|) is non-stationary away from the "wave front" |:r| = t. 
Hence 

sup \(G x ,ht)\<% t>l 

x:||x|/t-l|>e 1 
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for every integer n [IS]) an d it remains to estimate 



OO JJ. 

T ||x[i- £ ,i +£ ](k|/t)e-^^||. (24) 
o r 

Finiteness of (|24|) requires, in particular, that the electrons do not propagate at the speed 
of light, which is true in nature, but not precluded for the dynamics generated by the non- 
relativistic Schrbdinger operator H a f The easiest case occurs when the electrons are in a 
bound state ^ G RanE x (H), A < S, where, by ((H)) . 

supHe^le-^H < e^E x (H) < oo (25) 
t 

for some /3 > 0. Then obviously the integrand in ()24l) decays exponentially in time, and hence 
the limit (|21jl exists. If Vl/ is not in a bound state but its energy is insufficient for an electron 
to reach the speed 1 — e, then (|24j) is still finite, at least if H is the Hamiltonian ()37|) of the 
standard model |T7]. More precisely, (|2*4*)) is finite for all f in a dense subspace of RanE\(H) 
with A < E + m/2, and for e in (|24j) small enough. Here m/2 = mc 2 /2 is the non-relativistic 
kinetic energy of a particle at the speed of light. The assumption A < £ + m/2 thus ensures 
that no electron can reach the speed of light. 

The asymptotic field operators have the important property that 

a* + (g)RanE x (H) C RanE x+M (H) 

(26) 

a + (h)RanE x (H) C RanE X - m (H) 

if supp(g) C {k : \k\ < M} and supp(/i) C {k : \k\ > m}. Using ((2*6^1 and the existence of the 
limit (|21j) we prove existence of the limit ()2(Jj) and that 

lim e iHt a#(h ltt ) . . . a#{h n , t )e- iHt ^ = a* (ht) . . . a*{h n )^ (27) 

t— >oo 

if V 6 Ran£ A (#), A + Ej M j < s > where M j = sup{|A;| : ^(fc) / 0} and the sum £\ Mj 
extends over all creation operators in ()27|) . The main technical difficulty in this last step is the 
unboundedness of the asymptotic field operators (|21j) . 

From the above discussion it is clear that it is physically more sensible to describe the time 
evolution of the electron by a relativistic Hamiltonian such as 

H$ = V=K+1 + V(X), (28) 

in place of (^Q). Then (|24j) is finite for all ^ in any spectral subspace E\{H)7i with A < oo. 

ff 

Hence the asymptotic operators a+{h) exist on a dense subspace of H. The main results in 
Q71 Q21 US] apply to both electron-Hamiltonian, (JTJ) and 

To conclude this discussion of scattering states we remark that for Rayleigh scattering it 
suffices to prove existence of the asymptotic field-operators a+(h) on Ran-E/jx^g) (H ) , which 
follows from (|25|) . The improved results discussed thereafter are important in the study of 
photon scattering at a free electron (Compton scattering) |16| . 
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2.8 Asymptotic Completeness 

A characterization of ACR that is mathematically more convenient than Q19JI is achieved by 
mapping the freely propagating photons ha in (|18[) into an auxiliary Fock space T . We attach 
T to the Hilbert space TC by defining an extended Hilbert space TC = TC <8> T. The appropriate 
time evolution on TC is generated by the extended Hamiltonian H = H®l+l®Ht. Furthermore 
we define an identification operator J:DcH->Wona dense subspace D of TC by 

J* <g) | vac) = ^ 
J* ® a*(/ii) • • • a*(/i n )|vac) = a* (hi) • • • a*(h n )^ 

and linear extension. Since e~ -f*|vac) = |vac) and a*(ht) = e~ lH f t a*(h)e lH f t we can use I to 
write 

a*(hi, t ) • • • a*(h n , t )e- iHt V = Ie~ iflt [* a*(fci) • • • a*(/i n )|vac)] . 

If is an eigenvector of H and Pg denotes the orthogonal projection onto the closure of the 
span of all eigenvectors, it follows that 

a%(hi) ■ ■ ■ a* + (h n )^ = lim e iHt a*(h u ) ■ ■ ■ a* (h n , t )e- im V 
= n + [^f<S> a* (hi) ■ ■ ■ a*(/t n )|vac)] 

where 

n + = s - lim e iHt Ie- iHt P B 1 

t— +oo 

is the wave operator. Thus existence of scattering states becomes equivalent to existence of 
the wave operator Q + , and, since Ranf2 + = H+, asymptotic completeness as defined in (|19(1 
becomes 

%oo,S)WcRanl] + . (29) 

(It turns out that f2 + is a partial isometry and hence RanJ7 + is closed.) 

The reader familiar with quantum mechanical scattering theory is cautioned not to think of 
ACR as a form of asymptotic completeness for potential scattering. The comparison dynamics 
generated by H is not the free dynamics for all bosons. The bosons in the first factor of TC 
still fully interact with the electrons. Rayleigh scattering is more similar to iV-body quantum 
scattering with the additional complication that the number of particles is fluctuating. 

Asymptotic completeness of Rayleigh scattering, as described in Section 12.61 makes two 
assertions. First, any initial state ^ G 'RanE^_ 00 ^(H), that is not an eigenvector of H, in 
the course of time will relax to a bound state by emission of photons. Second, the asymptotic 
dynamics of the emitted radiation is well approximated by the free photon dynamics. Our 
proof of ACR contains two main technical ingredients that address these issues. In both of 
them we need to assume that either the photons are massive, i.e., u(k) = \/k 2 + m 2 , or that 
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an infrared cutoff a > is imposed on the interaction. In the second case m = a/2 in the 
following. 

Our first main ingredient is an estimate on the ballistic spacial expansion of the system for 
states with energy distribution away from S = a pp (H) + Nm. We show that S is closed and 
countable and that for each A G M\S there is an open interval A B A and a positive constant 
C\ such that 

(* t ,dT(y 2 )V t } >C x t 2 , t^oo (30) 

for all ^ C KanE/\(H). The proof is based on the positivity of the commutator obtained 
by differentiating the left hand side twice with respect to time. This positive commutator 
estimate, often called Mourre estimate, is proven by induction in energy steps of size m along 
a strategy very similar to the proof of the Mourre estimate for A-body Schrodinger operators 
|36| . We generalize the Mourre estimate in ^3] to accommodate our model. 

The second main ingredient is a propagation estimate for the asymptotic dynamics of 
escaping photons. Explicitly we show that 

f°° fit 

J j^ t JFdT(P t )Ff^ t )<C\\n 2 (31) 

for all \P E Ti, where 

P t = (Vlu - y/t) • X (\y\ > t S )(Vco - y/t) 

and < 5 < 1. Here / is an energy cutoff, F = F(&T{y 2 /t 2 A 2 )) a space cutoff, and A > a 
parameter that is chosen sufficiently large eventually. The left-hand side of (|31j) compares the 
average photon velocity, y/t, with the group velocity, Vw, for photons in the region {\y\ > t s }. 
This includes all photons that escape the electron ballistically. The finiteness of C thus confirms 
that the dynamics of outgoing radiation is approaching the free photon dynamics in the limit 
t — > oo. 

Asymptotic completeness had previously been established for a model with S = oo (con- 
fined electrons), and massive photons u(k) = \Jk 2 + m 2 , m > 0, by Derezinski and Gerard 
|14j . The methods in could probably be extended to prove ACR for our system. Instead 
of doing so, we chose to give an entirely new prove of AC based on the relatively elementary 
propagation estimate Q31|) . and using (|3U|) as the only dynamical consequence of the Mourre 
estimate. Our work is inspired by the Graf-Schenker proof of asymptotic completeness for 
A^-body quantum systems [22] . 

2.9 Absence of Excited States 

At present the knowledge on absence of eigenvalues above inf o~(H g ) is far less complete then, 
e.g., our knowledge regarding existence of a ground state. Known results on absence of eigen- 
values are derived under the assumption that g > is small enough [3 |H1 El , and to ensure 
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that no new eigenvalues emerge near inf &{H) an infrared cutoff is imposed 18 . There is a 
further assumption, the Fermi golden rule condition, which ensures that eigenvalues of H g= o 
dissolve for g ^ 0. This assumption can be checked in any explicitly given model. 

For the model introduced in Section 12. H with the assumption of an infrared cutoff, the 
following results hold true. For any given e > and for \g\ > small enough, depending on e, 

Vpp(H g ) n (mia(H g ), S at - e) = 0, 

where S at = infcr(-ff a t) jlUl I18j . This result, combined with ACR from the previous section 
implies that RanO + D HanEr^ jm^-^if) where the projector Pg in the definition of 0+ 
is the projector onto the ground state. That is, every ip G Ran.E( inf a (H),T, at -e)(H) relaxes to 
the ground state in the sense of Section 12.61 (|17jl . 

It has been asserted in [S| that the methods of can be used to show absolute continuity 

of the spectrum of H above and away from X a t for small \g\. This is presumably correct but 
a proof is missing. 

The strategy for proving absence of eigenvalues in a given spectral interval A C M is clear 
and simple: One tries to find a symmetric operator A on TC, such that 

E A (H)[iH,A]E A (H) > CE A (H) 

with a positive constant C. Since, formally, [iH,A]^f) = for every eigenvector of H, it 
immediately follows that a pp {H) n A = 0. The main problems, of course, are to find a suitable 
conjugate operator A, and to make these formal arguments rigorous. 

2.10 Relaxation to the Ground State 

An important consequence of AC for Rayleigh scattering and the absence of eigenvectors besides 
a unique ground state is relaxation to the ground state. To explain this let A denote the 
C*-algebra generated by all operators of the form 

B ® e^ h \ B G L(W at ), h G C^°(IR 3 ), 

where <fi(h) = a(h) + a*(h). We say that = e~ tHt ^ relaxes to the ground state ^o, if 

lim {* t ,A* t ) = (#o,^*o>(*,*> (32) 

for all A £ A. Suppose H has a unique ground state and let TL + denote the space of 
scattering states over ^o- That is, is the closure of the span of all vectors of the form 

a* + {hx)---a\(h n )^Q. 



16 



Then, by a short computation, all states in TC + relax to the ground state Vl/o |T7|. Since the 
assumptions of Section [2.91 imply AC in the form 

TL+ D Ran£^( inf(T ( H ) iSat _ £ )(i7 9 ) 

for given e > and small enough coupling \g\, it follows that all state in RanE^f -(H) i s at -e)(^g) 
relax to the ground state in the sense of Equation 



3 N-Electron Atoms and Molecules 

We now briefly describe how the results of the previous sections are generalized to the case 
o£ N > 1 electrons. For simplicity we neglect spin and Pauli principle. A (pure) state of N 
electrons is described by a vector ip £ L 2 (M? N ), and the Schrodinger operator for iV electrons 
in the field of K static nuclei is given by 



N 1 
H at , N = £ (-A Xj + V z { Xj )) + £ j— .. 



j=l i<j 

where Xj € M 3 is the position of the jth electron. The coupling of the electrons to the radiation 
field is done by a straightforward generalization of ( [o] ). The Hamiltonian of the entire system 
is given by 

Hg,N = #at,7V ® 1 + 1 ® H f + gH int 

N 

3=1 



9=0) 



and acts on L 2 (M. 3N ) (g) J 7 . Again, H g is self-adjoint on Z?(//, 
3.0.1 Spectrum and eigenfunctions of H at ,N 

Like .ff at , -ff a t,7V is a Schrodinger operator of the general form — A+V, and hence inf <7 ess (if a t i Ar) 
is given by Persson's theorem: 



inf a ess (H attN ) = lim inf ((p, H at N <p) i =: i, a t,JV 

R^co \tp£D R ,\\ip\\=l J 

where = Cq°(|X| > -R). Using the decay of the two-body potentials and the electron- 
electron repulsion one shows that X a t,jv = inf o"(i? a t,Ar-i), which leads to 

inf (T css (i? ati Ar) = inf cr(F at ,jv-i), (33) 

a special case of the more general HVZ-Theorem |481 I3b| . For Z > N — 1 the system described 
by -ff a t AT_i has a net positive charge and can bind at least one more electron. It follows, 
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by a simple variational argument, that inf a(H a t ; jsr) < inf <r(iifat,Ar-i) ) which, by (|33j) implies 
that inf <r(/f a t i Ar) is an eigenvalue of H at ^. In fact, ff a t,Af has infinitely many (discrete) 
eigenvalues below infcr ess 

(Hs± ; n) HH1- The continuous part of the spectrum of H at ,N is the 
interval [inf a eBS (H a t jv), oo) 3 and this interval may contain further eigenvalues below |53| . 
For a discussion of the structure of the continuous spectrum, the reader is referred to |36| . 

The results (jSJ), (|1U|) on the decay of eigenfunctions hold for Schrodinger operators in 
arbitrary dimensions, hence in particular for H a t jy. Better, non-isotropic exponential bounds 
are known too, but they are expressed in terms of a geodesic distance p(X) w.r.to a metric 
in M. 3N that depends on the spectra of the Hamiltonians -Hat.fcj for k < N (21 ESJ. Explicit 
expressions for p are known for N < 3, and for atoms under the (unproven) assumption that 
the ionization energy increases monotonically as the electrons, one by one, are removed from 
the atom j!2| . 

This concludes our discussion of -ff a t,JV an d we return to the composed system of ./V electrons 
and radiation. 

3.0.2 Exponential decay and ionization thresholds 

The ionization threshold E is the least energy that an atom or molecule can achieve in a state 
where one or more electrons have been moved "infinitely far away" from the nuclei. In an 
iV-particle configuration X G M 3Ar , one or more electrons are far away from the (static) nuclei 
if and only if \X\ is large. In this respect there is no difference between N = 1 and N > 1 
besides the dimension of the configuration space. Since this dimension is irrelevant for the 
proof of ()13|) . our result on exponential decay for N > 1 and its proof are straightforward 
generalizations of result and proof for N = 1. Let Dr := £ D(H) \ ^(X) = if |X| < R} 
and let 

Stv= lim f inf (^,H N ^)) . (34) 
fi-»oo \*e£>fl,IMI=i / 

Then for all real numbers A and j3 with A + (5 2 < E, 

|| (e m ®1)E X (H N )\\ < oo. (35) 

In the case of only one electron subject to an external potential V that vanishes at infinity, 
such as Vzi we saw that Etv=i = inf cr{H\ — V). The proper generalization to N > 1 is analog 
to the HVZ theorem for iV-particle Schrodinger operators. We show that 

S iV = mm{E%_ N , + E%,} (36) 

where E^, is the least energy of N' electrons with no nuclei present, Z = [21] ■ Like the HVZ 
theorem, (|33p . for the bottom of the essential spectrum of H at ,N, equation (|36|) requires the 
decay of the interaction between material particles with increasing spacial separation. While 
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this decay is obvious for the instantaneous Coulomb interaction, it is more tedious to quantify 
for the interaction mediated trough the quantized radiation field. The main problem in proving 
(J36|). however, is to control the error which arises when the field energy is split up into two 
parts, one associated with the N' electrons far out and one with the other N — N' electrons. 
This error is proportional to the number of photons, as measured by the number operator Nf, 
which in turn is not bounded with respect to the total energy and thus not under control. To 
deal with this problem we first proof (|36j) with an IR-cutoff a > in the interaction and then 
we show that Q36|l is obtained in the limit a —* |251 124j . 

The characterization (|36j) of the ionization threshold is important for proving that S^v > 
inf ct(Hn). 

3.0.3 Existence of a ground state 

The dimension of the electron configuration space is inessential for proving that (|16|) guarantees 
existence of a ground state. Therefore inf a{H^j) is an eigenvalue of whenever 

inf o-(Hjsf) < Sat. 

However, it is much harder to verify this condition for N > 1. The only easy case occurs 
for spatially confining external potentials where S^r = oo. In a tour de force Lieb and Loss 
recently showed that 

E N < min {E Z N _ N , + E%} 

for all atoms and molecules with Z > N— 1 |43j . Combined with (|36|) this proves that E?, < Sjv 
and hence that Efr is an eigenvalue of Hfj indeed. 

4 Concluding Remarks and Open Problems 

The results we have described on localization of the electron, existence of a ground state and 
existence of scattering states are established in [2U EI within the standard model of QED 
for non-relativistic electrons (see Appendix ^J. Asymptotic completeness is proved for the 
dipole approximation of that model, Hamiltonian (|59|) . an approximation that is physically 
reasonable for confined electrons |18j . We don't expect serious obstacles in proving ACR for 
the standard model (with IR cutoff), but to do so appears prohibitive in view of the additional 
work due to the interaction terms quadratic in creation and annihilation operators. The most 
important and most interesting open problem in connection with Rayleigh scattering is to 
prove completeness without IR cutoff. This has been done so far only for the explicitly soluble 
model of a harmonically bound electron coupled to radiation in dipole approximation and 
for perturbations thereof |51j. Steps toward ACR for more general electron Hamiltonians have 
been undertaken by Gerard [2U] . 
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The problem of the emitted low energy radiation can be understood as one aspect of the 
more general question of the intensity of the radiation in Rayleigh scattering as a function 
of the frequency. Experimentally, sharp spectral lines with frequencies u given by Bohr's 
condition are observed. This condition says that fko is the difference between the energies 
of two stationary states, that is, between two eigenvalues of H a f From quantum theory 
this phenomenon is expected to be a consequence of the smallness of a, which allows one to 
compute transition amplitudes in leading order perturbation theory. Rigorous work in this 
direction is currently being done by Bach, Frohlich and Pizzo @. One also expects that 
eigenvalues of H at show up as resonances in the spectrum of H g and that the eigenvectors 
of H g= Q are meta-stable states for the dynamics generated by H g , if g > 0, with a life-time 
inversely proportional to the resonance width. Both these expectations have been confirmed 
by work of Bach et al, and by Muck [13]. While the existence of resonances and meta- 
stable states is consistent with the experimentally observed spectral lines, it does not fully 
account for them. It remains to be shown that, for small a, first order transitions between 
meta-stable states dominate the process of relaxation to the ground state and hence that the 
intensity is largest for radiation obeying Bohr's frequency condition. A related question is the 
one about a confirmation of the correspondence principle within QED. By the correspondence 
principle, the frequency of radiation emitted by a highly excited atom agrees with the angular 
frequency of a classical point charge on the corresponding Bohr orbit. This principle together 
with Bohr's frequency condition determines the distribution of eigenvalues of highly excited 
states. A rigorous derivation of the correspondence principle would therefore confirm - but 
not prove - the domination of Bohr frequencies at least in the low energy spectrum. 

Many further questions arise once we allow for total energies above the ionization threshold 
S. Then the atom can become ionized and the dynamics of the removed electrons is close 
to the free one. The first task is thus to study the scattering of photons at a freely moving 
electron, the so-called Compton scattering. This has been done in jltjj . where we established 
asymptotic completeness for Compton scattering for energies below a threshold energy that 
limits the speed of the electron from above to one-third of the speed of light. To do so, we had 
to impose an infrared cutoff, for otherwise no dresses one-electron states exist. 

The natural next step is to combine Rayleigh with Compton scattering to obtain a complete 
classification of the long time asymptotics of matter coupled to radiation. This would include 
the photo effect as well as the occurrence of Bremsstrahlung. 

There are also very interesting and difficult open questions related to the binding energy 
£ — En even for N = 1. From Section [2.51 we know that 

E N =i < £ + infcr(-A + V). 
That is, if V — ► 0, the binding energy T, — E^ = i with coupling to the radiation field is at least as 
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large as the binding energy — inf a{— A + V) without radiation. Physical intuition tells us that 
this binding energy should actually increase due to the coupling to radiation: the radiation field 
accompanying the electron, by the energy-mass equivalence, adds to the inertia of the electron, 
that is, makes it heavier and thus easier to bind. This mass renormalization can explicitly 
be computed in the dipole approximation and this has been used to prove enhanced binding 
by Hiroshima and Spohn j^H]- Without dipole approximation the mass renormalization is not 
known explicitly and enhanced binding has been established so far only for small a |261 127j . 
It is an interesting and challenging problem to establish enhanced binding without dipole 
approximation and for arbitrary a and A. 

Once there are two or more electrons, one would like to know, first of all, whether two 
electrons attract or repel each another in our model of matter. Of course, equal charges repel 
each other but this argument neglects the effect of the quantized radiation field, which is 
attractive. Two charges close to each other will share part of their radiation field. Since this 
reduces the energy to produce it, binding is encouraged. The questions is thus whether this 
binding effect may overcome the Coulomb repulsion. 

A Non-relativistic QED of Atoms and Molecules 

The purpose of this appendix is to describe atoms and molecules within UV-regularized, non- 
relativistic quantum electrodynamics in Coulomb gauge. We shall also comment on suitable 
choices of units, on representations of the theory that avoid the use of polarization vectors, 
and on the dipole approximation. For further information the reader is referred to [71 I1H IT5] . 

A.l Formal Description of the Model 

To write down the model quickly and in a form familiar from physics books we shall be 
somewhat formal at first, using operator-valued distributions and avoiding domain questions. 

The Hilbert space of pure states of ./V electrons and an arbitrary number of transversal 
photons is the tensor product TL = Ti^t <£> T where 

Wat := A^L 2 ^ 3 ; C 2 ), T := ©£° =0 ® n s L 2 (R 3 ; C 2 ), 

<8>" =0 L 2 := C, and where ®"L 2 (M 3 ; C 2 ), n > 1, stands for the symmetrized tensor product of n 
copies of L 2 (M 3 ; C 2 ). The vector Q := (1, 0, 0, . . .) € T is called vacuum. The one-particle wave 
functions in H at and T are C 2 -valued to account for the two spin and the two polarization 
states of the electrons and transversal photons, respectively. 

The Hamiltonian of an atom or molecule with static nuclei is a self-adjoint operator in 7i 
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of the form 

N 1 

H = J2 Y~ & " ( ~* V ^ + V^A A ( Xj ))] 2 + aV R 1 + 1 ® fl>, (37) 
i=i m 

where <7j = (o~j iX , aj iV , aj, z ) denotes the triple of Pauli matrices acting on the spin degrees of 
freedom of the jth electron and Xj G M 3 is the position of the jth electron. The constant 
m > is the (bare) mass of an electron and a = e 2 /(he) = e 2 is the fine structure constant. 
In our units h = 1 = c. Another common form of H is obtained by using that 

[a ■ (-iV x + VaA A (aO)] 2 = (~iV x + y/aA A (x)) 2 + y/Za ■ B(x) (38) 

where B(x) = curL4(x). 

The operator Vr acts by multiplication with the electrostatic potential 

W = -E£j^ + Er^b < 39 > 

J=l 1=1 J 1<J J ' 

of the electrons in the field of K static nuclei with positions R\ , . . . , Rk G K 3 and atomic 
numbers Zi, . . . , Zk G Z+. We use the short-hands x = (x%, . . . , xn) and R = (R\, . . . , -Rat). 

The operators Hf and A A (x), for fixed i£l 3 , are operators on Fock space. Ht has been 
defined in Section f2.ll and A A (x) can be expressed in the form 

d 3 fc 
!\k\<A a/2| 

The polarization vectors B\(k) G C 3 , A G {1,2}, are orthogonal to the wave vector re and 
normalized 

e* x (k) ■ e^k) = 6^, s x (k)-k = 0. (41) 

In addition we assume that £\(tk) = £\(k) for all t > 0. The operators at (re) and a A (re) are 
creation- and annihilation operators in JF. These are operator-valued distributions, formally 
defined by a\(k)Q = for all k G M 3 , A G {1, 2}, and by the canonical commutation relation 



A A (x) = (2vr)- 3 / 2 £ / {e A (re)V^a A (re) + e A (re) e -^a A (fc)} . (40) 



A=l,2 



[a k (k),a* tt (q)] = S Xli S(k-q), [a*(k), a*(q)\ = 0. (42) 

A rigorous definition of Aa(^) will be given in the next section. The constant A > in (|40|) is 
the ultraviolet cutoff. Photons with \k\ > A do not interact with the electrons under the dynam- 
ics generated by H. This is nonphysical but necessary to define A A ^(x) on a dense subspace 
of T. For A = oo not even the vacuum would be in the domain of A(x). In fact, by a formal 
computation using the properties of a A (re) and a* x (k), \\Ai(x)Q\\ 2 = const Jj fc , <A Irel"- 1 ^ 3 re — > oo 
as A — > oo. 

In the QED of Feynman, Schwinger and Tomanaga, removing the UV cutoff requires a 
renormalization of mass, charge and field strength, a procedure that is mathematically not 
sufficiently well understood yet. 
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A. 2 Atomic Units and Perturbation Theory 

To work in the small-a regime it is convenient to choose the UV-cutoff A and the nuclear 
positions Ri € M 3 fixed on scales of energy and length where the units are proportional to 
the Rydberg energy ma 2 /2 = mc 2 a 2 /2 and the Bohr radius l/(ma) = h 2 /me 2 . We shall 
therefore rewrite the Hamiltonian in these units. It is instructive to begin by first scaling 
electron position and photon momentum independently. Let U : TC — > 7i be defined by 
(Ucp) n (x, fci,... , k n ) = rf , / 2 v? n l 2 (p n (rix,iik ll . . .,fj,k n ). Then 

N 1 - 2 

^- 1 uhu* =^2^—r- fa ■ (-»v XJ + ^/i4/ P (^i))] 2 

j=i lmT] ^ (43) 

OL 

+ —V R i„®l + l®Hf, 
rj/j, ' ' 

which is most easily verified using the definition of A^{x) given in the next section. In order 
that 2mrj 2 ^i = 1 and tj/jl = a we choose r\ = (2ma)~ 1 and /x = 2ma 2 . Next we express the UV 
cutoff and the nuclear positions in these units, that is we replace 

A//i -> A, R/t) -» R, (44) 

a non-unitary change of the Hamiltonian! Thus in the new units the Hamiltonian reads 

r i 2 

\a l ■ (-iVi + a 3/2 ^ A («^))J + V R ® 1 + 1 ® fT/ (45) 



AT 

E 

i=i 



where the dependence on a is concentrated in electron-photon interaction a 3j/2 A\(ax). The pa- 
pers by Bach et al. concern the Hamiltonian (|45|). many others concern (|37|) . When comparing 
results that are valid for small a only, one must keep in mind that these Hamiltonians are not 
equivalent, not even for atoms: the substitution A — ► Aa 2 , which occurs in (|44jl. corresponds 
to the change m i— » m/a 2 of the electron mass, as follows from (|43[) with r/ = = a 2 . 

A. 3 Fock-Spaces, Creation- and Annihilation Operators 

We next give a rigorous definition of the quantized vector potential A\(x) and we shall comment 
on the self-adjointness of H. In order to prepare the ground for the next section we define 
Fock space, creation- and annihilation operators in larger generality then needed here. A good 
reference for this section is 

Given a complex Hilbert space t) the bosonic Fock space over t), 

T = F(t)) = © n > 5 n (® n f)) (46) 

is the space of sequences ip = ((p n )n>o, with cpo S C, ip n S 5 n ((g> n f)), and J2 n >o ll^n|| 2 < oo- 
Here S n denotes the orthogonal projection onto the subspace of symmetrized tensor products 
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of n vectors in f). The inner product in T is denned by 

n>0 

where (p n ,'>Pn) denotes the inner product in <g) n [). We use J-"^ to denote the dense subspace 
of vectors ip G T with ip n = for all but finitely many n G N. 

Given /i € f) the creation operator a*(h) : .Ffi n C J 7 — > J 7 is defined by 

[a*(/i)(^] n = y/nS n (h ® </? n _i) (47) 

and the annihilation operator a(/i) : ^fl n C J 7 — > .T 7 is the restricted to J-gn of the adjoint of 
a*(h). The operators a(/t) and a*(h) satisfy the canonical commutation relations (CCR) 

[a(g),a*(h)] = (g,h), [a*(g), a#(h)} = 0. 

In particular, [a(h),a*(h)] = \\h\\ 2 which implies that ||a(/i)</?|| + \\tp\\ and ||a*(/i)y|| + \\ip\\ are 
equivalent norms. It follows that the closures of a{h) and a*(h) have the same domain. On 
this domain a*(h) is the adjoint of a(h) [111 Theorem 5.2.12]. The operator 

<j>(h) = -^(a(h) + a*(h)) (48) 

is essentially self-adjoint on It is useful to note that 

[4>{g),4>{h)]=i^m{g,h). 



In the case of QED, f) = L 2 (M 3 ;C 2 ) with inner product (g,h) = J2x=i 2 / 9\(k)h\{k) d 3 k 
and Ai(x) = 4>{G x> i), G x ,i G f), i = 1,2,3, being the components of 

G x (k,X) = ^e x (k)e- ik ^, (49) 

where n(k) = (2ir)~ 3 ^ 2 x\k\<A(k)- More generally we may allow k to be any real-valued, spher- 
ically symmetric function with n/y/uj G L 2 (R 3 ). In particular \k(— k)\ = which implies 
that [Ai(x),Aj(y)] = for all x,y G M 3 and i,j G {1,2,3}. 

The Hamiltonian (|37|) is defined and symmetric on the dense subspace 

V = [ A? =1 C °°(R 3 ; C 2 )] ^(Co 00 ^ 3 ; C 2 )) (50) 

where ^ fin (C^°(]R 3 ; C 2 )) is the space of vectors <p = (p n )n>o G .F fin with ip n = ® n Cg°(R 3 ;C 2 ). 
Let Hq = —A 1 + 10 Hf. Then Hq is essentially self-adjoint on P, self-adjoint on D(Hq) = 
D(—A (g) 1) fl D{1 Hf), and i7 — i?o is bounded relative to Hq. It follows that the closure of 
H\T> is defined on D{Hq) and symmetric on this domain. Since H is self-adjoint on D(Hq), 
according to Hiroshima |28| . we conclude that H is essentially self-adjoint on D. Alternatively, 
the Hamiltonian (|37|) may be self-adjointly realized in terms of the Friedrichs' extension of 
H\T>, since this operator is bounded from below, or, by using the theorem of Kato-Rellich for 
A/a small enough 9 . (aA small enough for the Hamiltonian (US)) .) 



24 



A. 4 Avoiding Polarization Vectors 

The fact that the polarization vectors are necessarily discontinuous as functions of k = k/\k\ S 
S 2 , by a well-known result of H. Hopf, may lead to annoying technical problems To show 
how these problems can be avoided we construct a representation of H that does not depend 
on a choice of polarization vectors ^3 ESI- This representation is based on a description of 
single-photon states by vectors in 

f) T := {h € L 2 {R 3 ;C 3 )\h{k) ■ k = for all k}, (51) 

the space of transversal photons. Let u : L 2 (M 3 ;C 2 ) — ► 1)t be the unitary map 

A=l,2 

where {£a}a=i,2 are the polarization vectors employed in the definition of H, an let U : T —* 
T(^)t) be defined by 

Ua*(h)U* = a*(uh), Ufl = Q. 
It follows that UAi(x)U* = (f>(uG Xji ) where 

(uG xd )(k) = ^ 4®e x (k) jex {kr 



11/2 



e 



(ej — • e^)), (53) 



{ei, 62, 63} being the canonical basis of M 3 . Note that <j){uG x j) is one operator and not a triple 
of operators even though k \— > uG x j(k) is a vector- valued function. 

The Hamilton operator = UHU* is the desired new representation of -ff. It has the 
form of H in ()37|) with the only difference that the form- factor of A(x) is now given by ()53|) . 
a function in C°°(1R 3 \{0}). 

By choosing other unitary mappings u from L 2 (1R 3 ;C 2 ) onto f)T one may define more 
equivalent representations of QED. For example the map 

u 2 : (hi,h 2 ) ^ Yl A k (54) 

A=l,2 

leads to the representation of H where the quantized vector potential is defined in terms of 
the form factor 

(u 2 G x>l )(k) = ^e-^( ei Ak), (55) 

the choice preferred in |42j . 

For the mathematical analysis of systems of electrons interacting with photons it is often 
necessary to localized the photons in their position space. That is, the photon wave function 
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h e t)T is mapped to J(iVk)h where J G Cq° (M 3 ). Now J(iVk)h [)t unless J = or 
/i = 0, and projecting J(iVk)h back to would destroy the localization accomplished by the 
operator J(zVfc). The solution to this problem is to work on the enlarged one-boson Hilbert 
space h ext = L 2 (M 3 ;C 3 ) = \\t @ fyz which also includes the space of longitudinal photons 
\]L = {h\k A h(k) = 0}. The Hilbert space for the entire system becomes W ex t = W a t <8> Ftf)) — 
TL F^l) and we define a Hamiltonian on TL cx t by 

# ext = H T ® 1 + 1 ® # /)L 
- 1 

= E a^[°-i • H v ^ + » 1/2 A(x))] 2 + «Vk + A> 
j=l 

where A/(x) = <f)(uG x A as above, but now uG x j is considered as an element of h ext = 
L 2 (IR 3 ;C 3 ). The fake longitudinal bosons from do not interact with the electrons and 
hence do not affect the dynamical properties of the system. However, by definition of H ext 
they contribute additively to the total energy and need to be projected out at the end of any 
analysis of the energy spectrum. 

To conclude this section we return to a more formal representation of A{x) by expanding 
photon wave functions in terms of ^-distributions. Let 5k{q) = S(q — k). We define, formally, 

af(k) = a*( ej 5 k ), a#(k) = {a* (k) , a* (k) , a* (k)) . 

From the expansion h(k) = Y^=i I f l j(l) e j^k(q) d 3 q and the (semi)-linearity of aft{h) we obtain 

3 

' ' d :i k 



a(h) = E / hj(k)aj(k)d 3 k = j h(k) ■ a{k) 
{h) = ^J h j (k)a* j (k)d 3 k= J h(k) • a*(k) d 3 k. 



3=1 
3 

a*(h) = ' 

3=1 

In particular, in the representation defined by (|5Hj) . 

A(x) = J 0^P(k){e lk - x a(k) + e- ik - x a*(k)} d 3 k 

where P(k) denotes the orthogonal projection onto the plane perpendicular to k. If (|55|) is 
used then 

A{x)= J 0^kA{e ik - x a(k) + e- ik - x a*(k)}d 3 k. 
A. 5 The Dipole Approximation 

In the dipole approximation of QED the quantized vector potential A{x) in the Hamilton 1)3 7 Jl 
is replaced by A(0). By 1)38)) the Hamiltonian (|3*7|) then reduces to 

N I 

H dip = ]T ^H V -, + ^ 1/2 A(0)) 2 + aV R + H f (56) 

3=1 
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where the interaction with the electron spin has dropped out. Without loss of generality we 
may now describe the electrons by vectors in the smaller space 7i a t = A^ 1 L 2 (R 3 ) of spin- less 
iV-fermion systems. 

The "constant" vector potential in (|56|) may be gauged away with the help of the operator- 
valued gauge transformation 

U = exp (^a 1/2 x i • ' ( 57 ) 

also known as Pauli-Fierz transformation. Since U(—i\7 x .)U* = — iS7 x — a 1 / 2 A(0), UA(0)U* = 
A(0), and 

N n 
UH f U* = H f + Va^xj -E(0) + a||K|| 2 (^\ 



where E(0) = —i[Hf,A(0)] is the quantized electric field, we arrive at 



N ( 1 _ \ 

VH^V* = £ ( - 2^ + ^ • ^(°) ) 

.7 = 1 



+ aV/? + iT/ + 1 1 as: 1 1 2 ^ ay J . 



(58) 



The dipole approximation seems justified when all electrons are localized in a small neighbor- 
hood of the origin x = 0, that is, when the total energy is below the ionization threshold. It 
then seems equally justified to drop the last term in (|58|) and to multiply x ■ E(0) with a space 
cutoff g £ Co°(lR 3 ); the later serves to ensure that the Hamiltonian H remains semi-bounded 
(after dropping the last term). This leads us to 

N / I _ \ 

^di P = Yl ( "^ A ^ + ^9{xj)x 3 ■ E(0) + aV R + H f , (59) 
j=i 

which is also called dipole approximation of 1)3 7|) . It has the advantage, over (|56|) . to be linear 
in creation and annihilation operators, which may simplify the analysis. 

The Pauli-Fierz transformation ()57|) is very useful in the analysis of the original Hamiltonian 
(fSTI) as well. Its effect is to replace A{x) by A(x) — A(0) = (p(G x — Gq) where 



\G x (k) - G (k)\ 



K (k) fjk-x 



^K-x _ y 



< |fc| 1/2 | K (A:)||x|. 



Thus the IR-singularity of the form- factor in UHU* is reduced by one power of \k\ at the 
expense of the unbounded factor \x\. This factor, however, is compensated by the exponential 
decay whenever the total energy is below the ionization threshold (see Section \2. 4(1 . 
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